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Abstract 

Pulsed laser deposition (PLD) is a popular growth method, which 
has been successfully used for fabricating thin films. Compared to con- 
tinuous deposition (like molecular beam epitaxy) the pulse intensity 
can be used as an additional parameter for tuning the growth behav- 
ior, so that under certain circumstances PLD improves layer-by-layer 
growth. We present kinetic Monte-Carlo simulations for PLD in the 
submonolayer regime and give a description of the island distance 
versus intensity. Furthermore we discuss a theory for second layer nu- 
cleation and the impact of Ehrlich-Schwoebel barriers on the growth 
behavior. We find an exact analytical expression for the probability of 
second layer nucleation during one pulse for high Ehrlich-Schwoebel 
barriers. 

Keywords: pulsed laser deposition, submonolayer growth, Ehrlich-Schwoebel 
barriers 



1. INTRODUCTION 

Pulsed Laser Deposition (PLD) is a growth method increasingly used for the 
fabrication of thin films fChrisey and Hubler (19841 )1. It is especially suited 



for the growth of complex multicomponent thin films, e.g. high temperature 
superconductors ||Cheung et al. (1993Q 1, biomaterials ||Cotell et al. (199^ )1, 



or ferroelectric films [ [Ramesh et al. (199 1| )]. A great advantage of PLD is 



the conservation of the stoichiometry of virtually any target material in the 
deposition. 
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Figure 1: (a) A simulated MBE-grown surface with D/F = 10 8 . (b) A 
simulated PLD-erown surface with D/F = 10 8 and / = 0.01ML. The 
coverage of both surfaces is 0.05ML. 



The main feature of PLD is that the target material is ablated by a 
pulsed laser and then deposited onto the substrate. Thus, in one pulse many 
particles arrive at the surface simultaneously The time between two pulses 
is of the order of seconds and the pulse length usually is of the order of 
nanoseconds [|Cheung and Sankur (19881) 1. Therefore we are going to neglect 
the pulse length. 

In this paper we present computer simulations showing several features 
of the growth morphology in PLD. In the next section a short description 
of the typical island distance for PLD in the submonolayer regime will be 
given. Thereafter the growth behavior of PLD under the influence of Ehrlich- 
Schwoebel barriers will be discussed. The aim of this work is to give an 
insight into possible reasons why in some experimental situations PLD has 
been shown to produce better layer-by-layer growth than ordinary molecular 
beam epitaxy (MBE) ||Jcnnichcs et al. (1996 )]. 



2. ISLAND DISTANCE 

The control parameters of PLD are the intensity J, which is the number of 
particles deposited in one pulse per unit area, and the diffusion-to-deposition 
ratio D/F. The average deposition rate is given by F = I /At, where At is the 
time interval between two pulses. The intensity is measured in monolayers 
(ML), and D/F is dimensionless, as the lattice constant is set to unity. 
The surface morphology depends sensitively on the pulse intensity which 
qualitatively can be seen in fig. |l[ One notices that in MBE there are few 
large islands relatively far apart whereas for PLD with a high intensity the 
surface is covered with many small islands. However, if one reduces the 
intensity to one particle per pulse, one would expect PLD to produce the 
same island morphology as MBE. This is indeed the case, although the two 
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situations are not exactly the same, as in PLD the deposition takes place at 
fixed times whereas in MBE it is probabilistic and therefore the time intervalls 
between two depositions have a Poissionan distribution. This difference, 
however, only influences the island morphology at high D/F, where finite-size 
effects are setting in [[Hinnemann (200"0| )1. In the simulations presented here 



only D/F below the finite-size region are considered. The simulations have 
been performed on a 400 x 400 square lattice and the island distance has been 
measured at 0.2ML coverage, when the island density reaches its maximum 
but coalescence does not yet set in. In the following the dependence of the 
island distance on the ratio D/F and on the pulse intensity I is investigated. 
For small intensities we recover the well-known power law for the island 
distance in MBE 

id - (j ) y (1) 

with the exponent 7 depending on the dimension of the surface, the island 
dimension, and the critical nucleus £*, i.e. the smallest stable island contains 
atoms [(Btoyanov and Kashchiev (1981|) ; |Schroeder and Wolf ( 1995) ]. For 



a two-dimensional surface, compact islands and a critical nucleus of i* = 1, 
one obtains 7 = 1/6. One should note that here the islands are not compact 
but fractal, as can be seen in fig. [I]. This is due to the fact that edge diffusion 
is not considered in the simulations. The exponent 7 can be determined from 
simulations such, that one monitors the number of nucleation events in a layer 
(oc lp 2 ) as a function of D/F. The value obtained in the present simulations 
is 7 = 0.17 ± 0.01. For PLD with large intensities the island distance obeys 
a different power law, however: 

h oc r v . (2) 

In this regime the island distance is independent of the parameter D/F, since 
the adatoms do not make use of their diffusion probability, as they find an 
island and attach to it in a much shorter time as the time they are allowed 
to diffuse between two depositions. The two regimes described by ([]]) and 
(0) are separated by a crossover at a certain intensity, where the number of 
deposited atoms is of the same order of magnitude as the adatom density. As 
the average adatom density n in MBE scales as n oc (D/F)- 1+2y , the critical 
intensity has to show the same scaling behavior 

, D . -1+27 

h « D ■ (3) 



F, 

It follows from ([[]), (Q) and (|3]) that the island distance can be represented 
as a scaling law [[Westerhoff et al. (2000|) 1 



where / is a scaling function with the asymptotic behavior 

~ * \ I — — const, for y <C 1 

fiy) \~ y -» fory>l. (5) 
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Figure 2: The scaled island distance for PLD versus the scaled intensity. 



Since Id does not depend on D/F for high intensities, the factor (D/F) 1 in 
(J|) must be compensated by the independence of /: I c v ■ (D/F) 1 = const., 
which together with (H) leads to 

~ (6) 



1 - 2 7 

The data collapse of the simulation results according to (|4]) is shown in fig. 
^. The exponent v obtained from the slope in the double-logarithmic plot 
is v — 0.26 ± 0.01 and thus in agreement with the predicted value v = 
1/4. Furthermore fig. § shows that the island distance indeed obeys the 
proposed scaling law. Thus there is a critical intensity which divides the 
intensity parameter space in MBE-like and different behavior. Below the 
critical intensity, the island distance in PLD equals or is very similar to the 
one in MBE, whereas above the critical intensity it differs significantly. 



3. GROWTH WITH EHRLICH-SCHWOEBEL 
BARRIERS 

In most experimental situations of PLD and MBE there are barriers to in- 
terlayer transport present, i.e. an atom experiences an extra barrier in ad- 
dition to the diffusion barrier when hopping down from an island fEhrlich 



and Hudda (19661 ); [Bchwoebel and Shipsey (19661) ]. This barrier is termed 



Ehrlich-Schwoebel barrier (ES-barrier) and is of the order of ~ O.leV for 
metals | Ruggcrone et al. (1997|) 1. Intuitively the ES-barrier can be explained 



as follows: In order to jump down an island edge, the atom goes through a 
position with a lower coordination than at the saddle point of diffusion on 
an island. The ES-barrier Ees is illustrated in fig. 0. A very useful measure 
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Figure 3: (a) Illustration of the Ehrlich-Schwoebel effect. The atom at the 
edge(black) has to overcome an additional energy barrier Ees t° hop down, 
(b) A model surface at the onset of second-layer nucleation. The calculations 
are done in the limit D/F — > oo. 



for its effect is the Schwoebel length, which is defined as 

Ees 



lES " 6XP [k^f) (7) 

and will be used as the control parameter in the following. The Schwoebel 
length is dimensionless, as the lattice constant is set to unity. In general 
the Ehrlich-Schwoebel barrier impedes layer-by-layer growth, as it forces the 
atoms to stay on the islands for a longer time. This increases the proba- 
bility that two atoms on an island meet and nucleate before they can leave 
the island separately. A model surface with second layer nucleation in the 
submonolayer regime is shown in the right panel of fig. In the following 
a description of the growth with ES-barriers in the submonolayer regime is 
given. It is important to describe second layer nucleation as it contributes 
to surface roughening and therefore has to be reduced as much as possible in 
order to achieve layer-by-layer growth. This concept of formulating a second 
layer nucleation theory and predicting the growth mode for a given parame- 
ter set has been successfully applied to MBE [[Tcrsoff et al. (19941 )1, and can 
be generalized to PLD as follows. 

Whether one obtains layer-by- layer growth depends on the time during 
deposition of the first monolayer until nucleation in the second layer starts. 
If second layer nucleation does not start significantly before the first layer 
is completed, one gets layer-by-layer growth, but if it starts well before the 
completion of the first monolayer, the surface will roughen ||Villain (1991| )1- 



If one compares PLD to MBE under these circumstances, the high intensity, 
which produces small islands, seems advantageous, as adatoms on the islands 
can leave them more quickly, provided they are small enough. On the other 
hand a high intensity means, that many particles are deposited on the surface 
simultaneously, which increases the probability for second layer nucleation, 
as two particles may be deposited on the same island and meet before they 
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Figure 4: (a) The second layer nucleation probability p nuc for / = 0.01 
versus island size A for different Schwoebel-barriers. From bottom to top 
the Schwoebel lengths are l ES = 10°, 10\ 10 2 , 10 3 , 10 4 , 10 5 . (b) Plot of p nuc 
versus A for two different intensities and for high Schwoebel barriers Les — 
10 3 ,10 4 ,10 5 . 



can both leave the island separately. 

Now a description for the second layer nucleation probability p nuc in PLD 
will be given. Let this quantity be defined as the probability that on a given 
island of size A at least one nucleation occurs during one pulse. As it has 
been shown in a similar analysis for MBE, the probability of second layer 
nucleation can be used as a starting point of a second layer nucleation theory 
|[Krug et al. ( 2000]) ] . It will be calculated making the following assumption: 



If two or more atoms are deposited simultaneously on the same island, they 
meet and nucleate before they leave the island. This assumption is true for 
high enough Ehrlich-Schwoebel barriers. Moreover we assume that in this 
case the nucleation happens before the next pulse arrives, which requires that 
D/F is high enough. Then the probability of second layer nucleation equals 
the probability that two or more atoms are simultaneously deposited on an 
island of size A during a pulse of intensity /. Thus we obtain for D/F —>■ oo: 

Pnuc = l-exp(-IA)(l + IA), (8) 

where A is the island area. In order to verify (|8|), simulations have been 
performed, where the second layer nucleation probability was measured in 
dependence of the island size. The simulations were performed on a 400 x 400 
lattice for various intensities I and Ehrlich-Schwoebel barriers, whose size is 
indicated by the Schwoebel length ([?]). 

In fig. | to the left the results for p nuc are displayed for the intensity 
/ = 0.01 and for different Schwoebel-lengths. One can see that the measured 
curves approach the predicted curve for increasing Schwoebel-barriers and 
that the agreement is good for Ies — 10 3 , when the Schwoebel-barrier is high 
enough to fulfill the assumption that atoms deposited on an island meet and 
nucleate instead of leaving it. The agreement is not convincing for island 
areas below 10, as here the discrete nature of the island sites plays a role. 
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The measured curves also deviate from the predicted curve for island sizes 
above 40. These deviations are due to coalescence, which is not accounted 
for in the calculation of p nuc and which for I = 0.01 already starts for island 
sizes below 50. This explanation is supported by the diagram in fig. |]to the 
right. The measured curves and the predicted curve for / = 0.005 agree much 
better, as for a lower intensity coalescence starts at a larger characteristic 
island size. In general one can see that the second layer nucleation probability 
p nuc can be described adequately with (||) for island sizes smaller than the 
characteristic island size at the onset of coalescence. 



4. CONCLUSIONS 

The island distance depending on the intensity I and D/F has been de- 
scribed by a scaling law. For intensities larger than I c the typical island area 
scales like A oc I~ 2u . Inserting this into (§) leads to the conclusion that 
pulsed laser deposition in the limit D/F — > oo should give worse results than 
molecular beam epitaxy for layer-by-layer growth. This is clear, as molecular 
beam epitaxy leads to layer-by-layer growth in this limit even for high, but 
finite ES-barriers. We emphasize, however, that strain effects and transient 
mobility of freshly deposited atoms have not been taken into account in the 
present investigation. 

For finite D/F the situation is less clear. In order to further analyze this 
aspect, one needs a description of the time- dependent island area which will 
be published elsewhere [[Hinncmann et al. ( 2000| )1 . If one compares the growth 



mode of PLD and MBE, one can already infer, that the use of PLD can only 
be advantageous for materials with a relatively high Ehrlich-Schwoebel bar- 
rier, as there the small islands have the advantage that atoms deposited on 
them can leave them with a higher probability before they nucleate on them. 
For growth situations with a very small Schwoebel-barrier, the diminished 
island size is no advantage, on the contrary, as it enhances second-layer nu- 
cleation compared to thermal deposition. The conclusion is that it does not 
help to just increase the intensity in order to obtain better growth results 
with a high Schwoebel-barrier. At some point the increased second layer nu- 
cleation probability outweights the advantage of small islands | JHinnemann 
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